Comparison of microscopic models for disorder in bilayer graphene: Implications for 
the density of states and the optical conductivity 
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We study the effects of disorder on bilayer graphene using four different microscopic models 
and directly compare their results. We compute the self-energy, density of states, and optical 
conductivity in the presence of short-ranged scatterers and screened Coulomb impurities, using both 
the Born approximation and self-consistent Born approximation for the self-energy. We also include 
a finite interlayer potential asymmetry which generates a gap between the valence and conduction 
bands. We find that the qualitative behavior of the two scattering potentials are similar, but that 
the choice of approximation for the self-energy leads to important differences near the band edge 
in the gapped case. Finally, we describe how these differences manifest in the measurement of the 
band gap in optical and transport experimental techniques. 



I. INTRODUCTION 

The current interest in bilayer grapheneii^ is motivated 
largely by the possibility of its application in the design- 
ing of novel electronic devices, as well as the opportunity 
to investigate the fundamental physics of massive chi- 
ral electrons in a condensed matter setting. In particu- 
lar, the opportunity to open a dynamically tunable band 
gap by electronic gating was predicted theoretically^.^ 
and verified in optical experiments^^^ , and this presents 
an opportunity not currently available in traditional 
two-dimensional electron systems^. However, transport 
measurementsiSr— demonstrate that the full nature of 
the effects of disorder induced in the graphene by its 
environment^^ is not currently resolved. This is a cru- 
cial issue since it appears that environmental disorder 
is the main limitation on the favorable properties of bi- 
layer graphene devices. In particular, there is a wide 
discrepancy in the size of the band gap extracted from 
transport measurements (which has obvious negative im- 
pact on the application of bilayer graphene in switching 
devices) and significant variation in the sub-gap conduc- 
tivity as a function of temperature. At low temperature, 
variable range hopping via midgap states dominates and 
produces a relatively small value for the gap, while at 
higher temperature, thermally activated transport be- 
tween the band tails is the predominant mechanism-"^. 
It has been suggested^^^ that both charged impurity 
disorder and short-range defect scattering play a role in 
transport properties of bilayer graphene and the pud- 
dling of electrons due to charged impurity disorder was 
also shown to be crucial in the understanding of capaci- 
tance measurements of dual-gated bilayer graphene^Sr^. 
A very recent work^- shows that the percolation trans- 
port gap and the spectral band gap could be significantly 
different in bilayer graphene due to the strong poten- 
tial fluctuation and density inhomogeneity induced by 
electron-hole puddles arising from random charged im- 
purities in the environment. Therefore, disorder is a key 
effect in many different measurements and a thorough 
understanding of the microscopic origin and effects of 



impurities and other scatterers is highly important. 

Previously, theoretical study of the role of disorder in 
bilayer graphene has been undertaken^l^— in which sev- 
eral types of disorder (such as resonant scatterers, lattice 
vacancies, short-range scatterers, and generic scattering 
processes) were investigated within a variety of theoret- 
ical schemes (such as the coherent potential approxima- 
tion, the Born approximation and self-consistent Born 
approximation, as well as other numerical techniques). 
The effect of disorder on the density of states (DOS), 
dc and ac conductivities, and other physical quantities 
were studied. However, to date, there has not been 
a systematic comparison of the most physical scatter- 
ing mechanisms - charged impurity disorder and short- 
range scatterers - for different approximation schemes. 
This is the topic of the current article. In this pa- 
per, we present a comprehensive analysis of the effect 
of charged impurity disorder and short-range scatterers 
on the DOS and optical conductivity of bilayer graphene 
in the presence of a band gap created by electrostatic 
gating. We compare these two scattering mechanisms 
which model the electron-impurity interaction, and two 
approximations to the self-energy induced by this interac- 
tion in the perturbation theory. Specifically, we consider 
the screened electron-impurity Coulomb interaction (CI) 
and the short-ranged interaction (SRI) modeled as a i5- 
function potential in real space as the model disorder 
potentials for our study. This paper is the first anal- 
ysis which treats these two scattering mechanisms on 
an equal footing and makes direct comparisons between 
them, since previous works have treated only different 
forms of short-ranged electron-impurity interactions, and 
we extend our analysis of the CI presented previousljs^. 
For both of the these electron-impurity interaction po- 
tentials, we describe the Born approximation (BA) and 
self-consistent Born approximation (SCBA) and high- 
light the differences between these two levels of approxi- 
mation in the perturbation theory. We also consider both 
the gapped and gapless situations. It should be noted 
that some of the numerical techniques which were pre- 
viously applied to the electron-impurity interactio n^^i^^ 
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technically go beyond the BA and SCBA and their appli- 
cation to the gapped case would be desirable. However, 
we believe that particularly the SCBA will capture the 
relevant physics near the band edge. 

Although our main interest is a comprehensive theo- 
retical formal understanding of disorder effects on bilayer 
electronic properties using a number of different models, 
a strong motivation for our work is the large discrepancy 
between the experimentally extracted transport gap in 
bilayer graphene and the theoretically calculated band 
gap obtained from band structure calculations. We also 
want to understand the effect of disorder on the opti- 
cal conductivity of bilayer graphene since extracting an 
optical band gap from the optical conductivity data is 
non-trivial in the presence of strong disorder (although 
it is routinely done in the literature, often incorrectly in 
our view). Our work should have implications for the 
bilayer density of states (and hence the transport gap) 
and the optical gap extracted from optical measurements. 
We also mention a particular advantageous feature of our 
theory is that we use the accurate four band model for bi- 
layer graphene throughout this work and retain the con- 
tributions from all four bands, avoiding various simpli- 
fications used in the theoretical literature dealing with 
bilayer graphene. 

The structure of this paper is as follows. In Sec. lU we 
describe the self-energy of electrons in bilayer graphene 
using the short-range and Coulomb scattering potentials 
for the BA and SCBA, and in Sec. IIIII we apply the re- 
sulting Green's functions to compute the DOS. In Sec. lIVI 
we describe the optical conductivity and then in Sec. |V] 
we focus on the highly important issue of the band gap 
and how it is extracted from various experimental tech- 
niques and describe the role of disorder in that context. 
Finally, we summarize our results in Sec. I VII and high- 
light the most important conclusions of our work. Two 
appendices contain additional results: A scheme to use 
the computationally efficient SRI to approximate the CI 
for low carrier density, and the derivation of the optical 
conductivity. 

We now introduce the single-particle theory which un- 
derlies our considerations of disorder. Throughout this 
paper, we use the four band tight-binding model for bi- 
layer graphene but neglect next-nearest neighbor hop- 
pings between the two layers. We label the electron states 
with the band index c, v for the conduction and valence 
bands, and with the branch index I, s for the low-energy 
and split bands. The dispersion of the four bands is given 
by the equation^''^ 



i^kA = i^x\le + i + i + bx^i + e{u^ + 7?) (1) 

where A — {v\,bx} labels the band, v = +1(— I) for the 
conduction (valence) band, and h = -1-1 (— I) for the split 
(low-energy) branch [see Fig. IHa)], and C = Tiv-pk. The 
inter-layer coupling in the tight binding formalism is 71, 
and Up is the Fermi velocity of monolayer graphene^i^. 
This band structure is isotropic in the wave vector and 
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FIG. 1. (a) Sketch of the low-energy band structure of bi- 
layer graphene in the presence of an external electric field 
which creates a band gap at low energy. The bands are la- 
beled with the index c, v for conduction and valence bands, 
and I, s for low-energy and split bands respectively. Diagrams 
corresponding to (b) the Born approximation and (c) the self- 
consistent Born approximation for the self-energy, and the 
definition of the Green's function in terms of the self-energy, 
(d) The Thomas-Fermi screening wave vector qtf as a func- 
tion of the Fermi energy. The band edge in the u = 0.2eV 
case is at Ep = 0.083eV. (e) The screening wave vector in 
the RPA as a function of wave vector in the intrinsic case. 

symmetric between conduction and valence bands. In 
the presence of an external electric field perpendicular to 
the plane of the graphene, an asymmetry is introduced 
in the potential on the upper and lower layers which gen- 
erates a gap between the cl and vl bands at low energy. 
This gap is parameterized by the energy u, correspond- 
ing to the energy difference between the two bands at 
k — 0. However, the minimum of the band gap is found 

at a wave vector /cgap = 2K^\f^^f^' ^^'^ gwen by 

-^gap ~ "^^71 / + 7i- This is known as the "sombrero" 
shape. We do not include various next-nearest neighbor 
interlayer hops which have been considered elsewhere^i^. 
We justify this by noting that, for example, the energy 
of the inter-band optical transitions near the K point are 
only modified by a few percent of the interlayer coupling 
71 when these hops are included. In a similar way, while 
there will be some small quantitative difference in the 
DOS due to inclusion of these higher order terms, the 
position of the band edge is almost unchanged and the 
disorder is by far the dominant effect. Therefore the sim- 
ple tight-binding model we outline is sufficient to obtain 
an accurate description of the DOS and optical conduc- 
tivity in the presence of disorder. 
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The two-dimensional CI in momentum space is 



yC(q) 



(2) 



where d is the distance of the impurities from the 
graphene surface, k is the effective dielectric constant of 
the system, and qs is the screening wave vector, which 
can in principle be a function oi q. In the extrinsic case 
[i.e. when there is a finite density of electrons or holes 
in the graphene so that Ep 7^ 0), the screening is taken 
into account in the Thomas-Fermi approximation^^ so 
that the q dependence of the screening wave vector is 
constant and 



9tf 



-Dq{Ey) 



(3) 



where Do{Ef) is the non-disordered D OS at the Fermi 
energy. This screening wave vector is shown in Fig. Wid) 
as a function of the Fermi energy. In the ungapped case 
(dashed line), the DOS is finite for all values of Ep so 
qxF is well defined. The shape of the DOS is replicated 
in gxF so that it increases linearly except for the onset 
of the split bands which causes a step in the DOS at 
Ep « O.SeV. Note that this linear increase is a result 
of the more accurate hyperbolic model which we use for 
the dispersion of bilayer graphene near the K point. The 
quadratic approximation is only valid at low density and 
predicts a constant DOS, however the density ranges we 
consider exceed the applicability of this model-^. When 
the gap is present, the DOS is a more complex func- 
tion which is refiected in gxp [solid line in Fig. [ijd)]. In 
particular, it is divergent at the band edge so that the 
screening becomes very strong for low carrier density. For 
intrinsic bilayer graphene (that is, when the carrier den- 
sity is zero so that Ep — 0) in the gapped regime, the 
density of states is not defined and some other approx- 
imation to the screening wave vector must be made to 
avoid the unphysical divergence of the Coulomb interac- 
tion at g = 0. In Fig. [Ije), we show the screening wave 
vector as a function of q in the random phase approxi- 
mation (RPA) for intrinsic bilayer graphene. This shows 
that for u = Q the q — value is finite, which is consis- 
tent with an earlier analysis of the screening in bilayer 
graphene^ computed using a quadratic approximation 
for the band structure. For u — 0.2eV, the \ow-q behav- 
ior is almost linear with the q — value being zero. In 
both cases, the screening wave vector depends linearly 
on q at high wave vector. Therefore, the g — >■ limit 
shows that this more sophisticated approximation also 
fails to remove the divergence in ^'^(q) for gapped bi- 
layer graphene. In order to make an order-of-magnitude 
estimate for the screening in this case, we take the ex- 
trapolation of the higher-energy screening wave vector 
to g = [as shown by the dotted line in Fig. [Hd)]. For 
u = 0.2eV, this means gxpi ~ 0.5 where a — 0.246nm 
is the lattice constant. The precise value chosen will not 
significantly affect the DOS or optical conductivity. The 
regularization of the CI at g = is necessary for avoiding 



and artificial divergence, but the details of this infrared 
regularization do not affect the results presented in this 
paper. 

The short-range interaction (SRI) is given by the po- 
tential 



(4) 



where the interaction strength is parameterized by Vq. In 
principle, the parameter Vq is not known so it can either 
be used as a fitting parameter, or some physical reason- 
ing can be used to estimate its value. We make a phe- 
nomenological assumption that Vo — iire^dsc/ k where 
dsc is the scattering cross-section of the short-ranged im- 
purities. We assume that dsc = Inm and that the sample 
is mounted on an Si02 substrate so that k = 2.5 is the 
effective dielectric constant of the environment, giving 
Vb = 3.62eVnm2. In Appendix E] we demonstrate a way 
in which the SRI can be used to approximate the CI in 
the low-energy region. Our models of long- and zero- 
ranged disorder give similar results at low energies (see 
Appendix |^ so our model also applies in this energy 
range to the intermediate-range resonant scattering dis- 
order. In the literature, it is noted that the resonant scat- 
terers are likely to also contribute to the existence of mid- 
gap states in the system which, as far as we know, have 
never been observed in any optical experiments. We do 
not take these mid-gap states into account in our model, 
but mention that it should be straightforward to include 
such states in our theory if experiments demonstrate their 
existence. 



II. SELF-ENERGY 

In this section, we describe the self-energy introduced 
by interaction between the electrons and the impurities 
in bilayer graphene. We use two approximations, the BA 
and the SCBA, and compare the results of each. For- 
mally, the self-energy in the BA is 



|y(k-k')|2FkA,k'A' 



k'A' 



E-Ew 



k'A' 



■ ir] 



and in the SCBA is 



^SCBA 
^kA 



(E) = m 
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k'A' 



\V{k^k')\^F, 



kA,k'A' 



E-E, 



k'A' 



(5) 



(6) 



which differs from the BA because it includes the full 
Green's function on the right-hand side and therefore de- 
fines a self-consistent equation for the self-energy. The 
SCBA self-energy is generally found by iterating this 
expression taking the BA as the initial value of the 
self-energy on the right-hand side until convergence is 
reached. The diagrams corresponding to these approxi- 
mations are shown in Fig. [TJb,c). Although the SCBA is 
a self-consistent approximation involving an infinite se- 
ries of diagrams, it is not always necessarily better that 
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FIG. 2. The real and imaginary parts of the self-energy in the low-energy conduction band (solid lines) and split conduction 
band (dashed) hues for ni = lO^^cm"^ (red lines) and ni — lO^'^cm"^ (black hues) at |kta = 0.01. The scattering potential (CI 
or SRI), approxmation for the self-energy (BA or SCBA) and band gap (u = or u = 0.2eV) are labeled in each panel. The 
self-energy is trivially zero in the clean limit. 



the BA since in each order many diagrams (e.g. the 
crossing diagrams) are left out - in general, the self- 
consistency in the SCBA tends to smooth out the sharp 
features in the BA results as we shall see below in our 
calculations of the DOS. In these equations, the symbol 
rii is the areal density of charged impurities, V^(k) is the 
electron-impurity interaction potential for the wave vec- 
tor k, FkA,k'A' is the wave function overlap between kA 
and k'A' states, and is a positive infinitesimal. In the 
following we shall describe the self-energy of both the CI 
and the SRI each of which requires the appropriate ex- 
pression for the interaction potential [Eq. ^ or Eq. ^] 



to be substituted into the equations for the self-energy. 
Because the bands are electron-hole symmetric, we also 
have the symmetry relation 

Ekc.(i?) = ^Kvui~E)- (7) 

The real and imaginary parts of these self-energies are 
plotted in Fig [51 Specifically, we take wave vector 
|k|a = 0.01, impurity concentrations rii — lO^^cm^^ and 
rtj = lO^^cm"^, K — 2.5 corresponding to a back-gated 
Si02 sample, 77 — 10~'^eV and two values of the gap: 
u = and u — 0.2eV as indicated by the labels in Fig. [21 
The solid lines refer to the Ic band, the dashed lines to 
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the sc band. The Iv and sv bands can be found via 
Eq. ([7]). In the CI case, we have taken n = 5 x lO^^cm^^ 
to define the screening wave vector. Throughout this 
article, we have 71 ~ 0.30eV, 70 — 3eV which gives 
VP = y^-foa/i2h) = 9.75 x 10'5ms"\ and = 0. In the 
SRI/BA, the factor UiVq is a multipher, meaning that the 
fines for the different impurity concentrations are scaled 
copies of each other, but this is not the case for the CI or 
for the SCBA. We see that the BA and SCBA give rather 
similar results for most values of energy, especially in the 
gapless case. However, the sharp spikes in the BA due 
to band edges are rounded off in the SCBA, and the on- 
set of the imaginary part at the band edge is shifted to 
slightly lower energy in the SCBA. These differences will 
be manifest in the DOS at the band edge as we discuss 
below. In the u = case, the imaginary part of the low- 
energy branch is always finite, and is roughly constant 
up to the energy where the split bands become occupied. 
When u is finite, the imaginary part disappears in the re- 
gion near E = signifying that a gap has opened. Large 
peaks appear near the band edge, and these peaks are 
not symmetrical in energy, but are larger on the positive 
energy (conduction band) side of the gap. This behav- 
ior is reversed in the valence band due to the relation 
in Eq. ([7]). The presence of the gap also modifies the 
intcr-band contribution to the self-energy, reducing the 
size of the imaginary part in the higher energy regions. 
In addition to the change in the relative size of the self- 
energy, increasing impurity concentration has the effect 
of introducing a larger shift to the energy at which the 
band minima occur. 



III. DENSITY OF STATES 

These self-energies described in Section |ll] define 
Green's functions which can be used to compute the DOS 
as 

D{E) = -^-^y2lmGME) (8) 

TT ^ — ' 
kA 

where and are the spin and valley degeneracies 
respectively, and G^xiE) ^ [E - E^x - ^kxiE)]"^ is 
the Green's function incorporating the self-energy in- 
duced by the electron-impurity interaction (see Fig. [IJ . 
The DOS for the clean system is found by substitut- 
ing 'SkxiE) — —ir] (where 77 is a positive infinitesimal) 
or by exact extraction from the band structure. The 
DOS functions are shown in Fig. [3] for the BA and the 
SCBA, for the CI and SRI, at two different values of im- 
purity concentration, and for the ungapped and gapped 
cases. In the gapless situation (top row of Fig. [3]) the 
only difference between the two approximations is a slight 
rounding-off of the step at the onset of the split bands in 
the SCBA. For our choice of parameters, the SRI gives a 
stronger change in the slope of the DOS with increasing 
impurity density than the CI, leading to a larger change 



in the DOS at higher energies. Also, increasing impu- 
rity concentration shifts the bottom of the split bands to 
lower energy, as indicated by the position of the step-like 
feature. In the gapped case (middle row of Fig. in- 
creasing impurity concentration has the marked effect of 
lowering the energy of the onset of the low-energy band. 
In fact, if the disorder is sufficiently strong (and the clean 
band gap is not too large), the gap can be closed by dis- 
order broadening''''. Also, when the gap is present the 
features of the band edge are qualitatively different be- 
tween the BA and the SCBA. Specifically, the divergence 
present in the clean DOS at the band edge persists as a 
sharp spike in the BA, but is rounded off in the SCBA, in- 
dicating that the DOS approaches zero in the low-density 
limit for the SCBA but continues to diverge in the BA. 
Probes which are sensitive to the DOS, such as the tun- 
neling conductance in atomic force microscopy, capac- 
itance measurements, or direct measurement of ^ via 
single electron transistor spectroscopy should be sensitive 
to this difference at low carrier density. To highlight this 
difference, in the bottom row of Fig. |3]we show the low- 
energy conduction band edge for three different values of 
the impurity density for each pair of self-energy and scat- 
tering potential approximations. In the BA [Fig. [3Uc,i)], 
increasing disorder decreases the height of the peak at the 
band edge, and creates a dip in the region where the band 
edge would be in the non-disordered case which widens 
with increasing impurity density. The SCBA shows very 
different behavior at the band edge. The sharp peak 
disappears, and the onset of the finite density of states 
occurs at a slightly lower energy than in the BA. The CI 
and the SRI give similar results for the modification of 
the band edge (except that the SRI is stronger for equiv- 
alent choice of rii because of our choice of Vq) indicating 
that the SRI can even be used to reliably approximate the 
low-density behavior of the DOS in the CI case as long 
as Vo is set correctly, as described in Appendix [A] The 
reason for this is that the screening wave vector is of the 
order of the size of the Brillouin zone [see Fig[TId)] which 
is large in comparison with the wave vectors near the K 
point which dominate the contribution to the Green's 
function at low energy. Therefore, the q dependence in 
V'~' is negligable at this energy scale and the interaction 
strength is almost independent of q, just as it is in the 
short-range case. However, for all impurity densities and 
for both approximations, the DOS remains close to the 
non-disordered value away from the band edge with the 
disorder broadening affecting the DOS mainly near the 
band edge. 

IV. OPTICAL CONDUCTIVITY 

We now turn our attention to the optical conductivity 
of gapped bilayer graphene. This is an important quan- 
tity from the experimental perspective since the optical 
identification of graphene^ii^^ and the optical determi- 
nation of the band gap^""— both depend crucially on a 
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FIG. 3. The density of states for ni = lO^^cm"^ (red lines), rii = lO^'^cm"^ (black lines), and (near the band edge) for 
Tii = 5 X lO^^cm"^ (blue lines). The scattering potential (CI or SRI), approximation for the self-energy (BA or SCBA), and 
gap {u = OT u — 0.2eV) are labeled in each panel. The dotted lines give the DOS in the clean limit. 



thorough understanding of the optical conductivity. The 
optical conductivity in the absence of disorder was cal- 
culated some time ago^i^, but we extend this analysis 
to include the effects of disorder using the same CI and 
SRI models discussed before. The optical conductivity is 
given by 

.^x . f k'dk' /"^^ dE' 

X Mly{k')lmG^.x{E')lmG^.y{E' + E) (9) 

where 

MlAk)^j^ ^|(A,fc,O|n{;,|A',fc,0)^ (10) 

The derivation of these equations is given in Appendix 
IbI We use = ImeV and show the results in Fig. |4] 
for the four combinations of scattering mechanism and 
approximation for the self-energy. Figure |4l[a)-(d) show 



the intrinsic case, and the origin of each of the features 
in the curves is sketched in Fig. SJ^i). In the CI shown in 
Fig.mja-b), the peak associated with the interband tran- 
sitions is shifted to lower energy with increasing disorder. 
In the BA [Fig. Hl^a)] the sharp spike at the band edge 
persists even in the presence of strong disorder, while in 
the SCBA, [Fig. IHJb)] it is rounded off and a significant 
tail develops. Also, in the BA, there is a small peak near 
the energy E^^^ = 166meV which is an artifact of this 
approximation and is not present in the SCBA. These 
features directly reflect the related structure in the DOS. 
Comparing with the short-range interaction in Fig. Sfc- 
d), it is clear that the qualitative behavior of the optical 
conductivity near the band edge is identical to the CI. 
However, because the strength of the interaction remains 
constant at large wave vectors, the structure near the 
onset of the transitions including the split band is very 
strongly modified even for small values of the disorder. 

We have taken the carrier density of n = 5 x lO^^cm^^ 
as an example of the structure of the optical conductiv- 
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FIG. 4. Panels (a-h) The optical conductivity of bilayer graphene in units of ao — (the conductivity of monolayer graphene 
in the clean limit) with u — 200meV. The type of scatterer (CI or SRI), approximation of the self-energy (BA or SCBA), and 
carrier density are labeled in each panel. In all plots, the dotted line represents the contribution from the interband conductivity 
in the clean system (i.e. the Drude peak is not included in the extrinsic case), the black line is ni — lO^^cm"^, and the blue 
line is rii = 5 X lO^^cm"^. For the finite disorder cases, the intraband contribution is included. Throughout, we have assumed 
a Si02 substrate with no top gate when defining the dielectric environment (see text for full definition of parameters). Panels 
(i), (j) show the transitions which give rise to each feature in the non-disordered conductivity. 



ity in the extrinsic case. Figures Sl^e-h) show the results 
for the same parameters as for the calculations of intrin- 
sic BLG. At this density, the Fermi energy is well above 
the sombrero region and therefore the effects of this low 
energy structure are not present. The features of the in- 
terband part of the optical conductivity are sketched in 
Fig. [3I^j). The additional structure in comparison to the 
intrinsic case is due to the Fermi-blocking of some transi- 
tions, and the additional transitions from the low-energy 
conduction band. The large peak at low energy is the 
intraband contribution (called the Drude peak), and is 
absent in the intrinsic case. The most obvious conclu- 
sion is that all the structure that is due to the non-equal 
transition energies is very quickly made indistinguishable 
by disorder, and even at rij = lO^^cm"^, the structures 
associated with the different interband transitions which 
include the split band become blurred into a single peak. 
Finally, for strong enough disorder, the interband tran- 



sitions become blurred with the intraband peak making 
identification of the individual features of the conductiv- 
ity impossible. Note that as the carrier density varies, the 
energies associated with transitions 1 and 2a will change, 
strongly modifying the optical conductivity. Therefore, 
accurate control of the density is essential for extracting 
the correct optical conductivity. 

V. BAND GAPS 

We now describe the effect of disorder on the band gap 
as measured by optical experiments in bilayer graphene 
and compare to the gap predicted by the DOS calcu- 
lations. There have been several attempts to determine 
the band gap via measurement of the optical conductivity 
in absorption or reflection measurements^^—. Figure 
shows the transitions in the intrinsic case. The transition 



8 



0.5 



^0 

00 

0.5 




0.5 1 1.5 2 



(a) CI/SCBA 



Optical 
DOS 



(c) SRI 

-Optical, SCBA 

DOS, SCBA 

DOS, BA 

Optical, BA 




15 



5 10 
"imp (10^^ cm"^) 



Optical gap 




20 



DOS gap 




FIG. 5. The band gaps extracted from the optical conduc- 
tivity of intrinsic bilayer graphene and the DOS. (a) SCBA 
approximation for the CI on an Si02 substrate; (b) the same 
approximation but for a high quality hBN substrate (note the 
difference in the range of the horizontal axis); (c) SRI for BA 
and SCBA; (d), (e) definition of optical and DOS gaps. 



between the low-energy valence and conduction bands is 
labeled as transition 1. This transition gives a very clear 
peak in the optical conductivity, and since there is no 
Drude peak in the intrinsic case, it serves as a clear way to 
define the band gap. The onset of the peak occurs at E^p 
in the clean limit, and therefore serves as the definition 
of the gap in this case. This definition is shown schemat- 
ically in Fig. ^d). Figure |4l^j) shows the transitions in 
the extrinsic case, and we refer to this sketch in order 
to explain the procedure for extracting the gap. The di- 
rect transition between the low-energy conduction and 
valence bands [labeled transition 1 in Fig. HJ^j)] is prob- 
lematic for the extraction of the gap because the energy 
of the transition is Fermi-blocked for the lowest energy 
excitations and therefore does not directly give the band 
gap. Some experiments^ instead computed the gap as the 
difference in energy between transition 2b and transition 
3, but this is also unreliable since even a small amount of 
disorder masks the intricate structure of the optical con- 
ductivity. Other experiments^!^ fit single particle theory 
with phenomenological disorder broadening to the exper- 
imental data, and while this procedure is more reliable 
than identifying particular peaks with individual tran- 
sitions, it cannot take into account the specifics of the 
impurity-electron interaction, particularly at the band 



edge. For the purposes of the following discussion, we 
define the energy of the optical band gap as the top of 
the first peak in the optical conductivity in the intrinsic 
case, as indicated in Fig. \5ld). Figure [5l^ a) shows that 
this point is very robust against disorder for the CI in the 
SCBA [compare to the position of the top of the peak in 
Fig. IH^b)]. The peak in the BA moves to lower energy 
with disorder, implying that the measured optical gap in 
this approximation is not robust compared with SCBA. 
Similarly, the short-range interaction for the SCBA leads 
to such strong modification of the peak structure that 
it is difficult to define the optical band gap at signifi- 
cant impurity concentration. This information is shown 
in Fig. [SJa) where the optical band gap is plotted as a 
function of impurity density. In the main panel, a Si02 
substrate is assumed. The solid lines denote the optical 
gap, and we see that for the CI in the SCBA [solid black 
line in Fig. [SJa)], the gap docs not close. In contrast, 
for the SRI [solid black line in Fig. ^h)] it is only possi- 
ble to define the band gap for a small range of impurity 
concentration, and the SRI in the BA [solid gray line in 
Fig. [SJb)] decreases rapidly. 

We compare this to the DOS gap, which we define as 
the energy difference between the onset of the DOS asso- 
ciated with the conduction and valence bands as sketched 
in Fig. El^e). This is plotted in Fig. [5Ua,b) using the 
dashed lines as labeled in each plot. It is clear that the 
effect of disorder is to rapidly shrink the size of the DOS 
gap in all approximations. For reference, we point out 
that using very high quality substrates, such as hexagonal 
boron nitride (h-BN)iii^, where the impurity concentra- 
tion is up to an order of magnitude smaller than in Si02 
substrates, the difference between the optical and DOS 
gaps defined above will be significantly smaller. This is 
demonstrated in the inset to Fig. ^a) , where the optical 
and DOS gaps for bilayer graphene on h-BN are plotted 
for the CI in the SCBA, for a reduced range of impurity 
density. 



VI. SUMMARY AND CONCLUSION 

In summary, we have investigated the effect of scatter- 
ing from charged impurities and short-range potentials 
on the self-energy, density of states, and optical conduc- 
tivity in gapped bilayer graphene using both the Born 
approximation and the self-consistent Born approxima- 
tion. We have demonstrated that these two approxima- 
tions give qualitatively different results for the DOS at 
the band edge, where the SCBA predicts that the di- 
vergence of the DOS seen in the clean limit and in the 
BA will be rounded off. Both approximations for the 
self-energy yield a DOS where the gap between the con- 
duction and valence bands closes with increasing disor- 
der strength. The main differences between the CI and 
the SRI occur at the onset of the split bands. For our 
choices of parameters, the SRIs are stronger at these en- 
ergies and therefore the shift in the band is larger than 
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for CI. We also introduced a self-consistent procedure to 
provide a computationally advantageous approximation 
for the low-energy properties in the CI via the SRI. We 
then demonstrated that the optical conductivity in the 
presence of disorder is a highly complicated function, the 
qualitative features of which depend crucially on the dis- 
order. Our two main conclusions are that even a small 
level of disorder is sufficient to blur the main features 
of the optical conductivity of extrinsic graphene, and 
this makes the determination of the band gap via this 
measurement highly problematic. In contrast, the peak 
associated with the interband transition of the optical 
conductivity in the intrinsic case is rather robust against 
the effects of disorder in the CI/SCBA case. Finally, we 
discussed the evolution of the band gap in the DOS and 
extracted from optical measurements of intrinsic bilayer 
graphene as a function of impurity density, finding that 
the DOS gap is quickly reduced with increasing disor- 
der, and eventually will be closed completely. In con- 
trast, the gap extracted from the optical conductivity of 
intrinsic bilayer graphene in the SCBA is much more ro- 
bust against disorder. This finding, that the DOS gap 
is always substantially smaller than the nominal optical 
gap (see Fig. [5]) defined operationally for gapped BLG in 
the presence of disorder, may provide an explanation for 
why the experimental transport gap extracted for gapped 
BLG appears to be substantially smaller than both the 
theoretical single particle band gap and the experimen- 
tally extracted optical gap^"— i^"— . Our results clearly 
indicate that this discrepancy between transport and op- 
tical gaps should disappear in systems with very little 
disorder such as high mobility suspended graphene or 
high quality graphene on h-BN substrates. It is also in- 
teresting to note that our conclusion regaring the large 
difference between transport and single particle (or op- 
tical) gaps is completely consistent with a recent theo- 
retical conclusion based on percolation considerations of 
BLG graphene in the presence of electron-hole puddles^. 
We thank US-ONR and NRLSWAN for support. 



Appendix A: Approximation of Coulomb scattering 
by short-range scattering 

In this appendix, we discuss a method for approximat- 
ing the low-energy properties of the CI using the SRI. 
The physical mechanisms of the SRI and the CI are not 
related to each other since different types of disorder in 
the environment of the graphene will lead to each. How- 
ever, as we have shown in this paper, the SRI and CI give 
qualitatively similar behavior for the physical observables 
such as the optical conductivity. If one is interested in 
transport characteristics, the compressibility of the elec- 
tron liquid, or other quantities that depend on the carrier 
density, the CI is a much more time-consuming approxi- 
mation to use because the screening properties vary with 
the carrier density n and hence the self-energies also de- 
pend on it. This implies that the SCBA procedure must 
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direct comparison of the DOS for CI and SRI for u = 0.2eV. 
The impurity density is labeled in each plot. In (c), V' = 
2eVnm^ and in (d) V = l.BeVnm^ 



be implemented for every value of n. This is not the 
case for the SRI, where the impurity interaction is con- 
stant with n and so one self-energy function is valid for 
all carrier densities. Therefore, it would be advantageous 
if we could describe a method for approximating the CI 
by suitable choices of parameters for the SRI. In this ap- 
pendix, we denote the value of Vq set with reference to 
the CI by the notation V'. The obvious method for doing 
this is to set 



V' = V^{0) 



which is still a function of n via the screening wave vector. 
If this function is reasonably flat then a suitable constant 
value can be chosen as an approximation and this step 
removes the dependence on n from the self-energy. The 
dashed lines in Fig. [5] (labeled 'non-renormalized' in the 
legend) show this for u = and u = 0.2cV. We see that 
for higher carrier density (n > 2 x lO^^cm"^) there are 
no sharp changes in the value of V' determined by this 
method, and that the slope is relatively modest. How- 
ever, at low carrier density in the gapped case, the inter- 
action strength becomes small and this approximation 
will not work. 

A further sophistication can be introduced to account 
for the effects of disorder on the density of states and 
hence on Qs- In this approach, the disordered DOS is 
computed in the BA using the value of V' obtained above, 
and the renormalized Fermi energy Ep is extracted by 
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solving 



E' 



D[E)dE = n 



where is a time-ordering operator for an imaginary 
time r, Vn = is a Matsubara frequency and (3 = 
l/ikBT). 

Then the lowest-order correlation function is given by 



where we assume that the disorder does not change the 
carrier density. This Fermi energy and DOS function can 
be used to set a new value of V' , and this procedure is 
iterated until convergence is reached for V'. Doing this 
for two disorder concentrations in the gapped and gap- 
less cases are given in Figs.[n^a,b). These plots show that 
selecting V' ~ 2eVnm^ will be a reasonable approxima- 
tion for higher carrier density. In Figs. [6][c,d) we show a 
comparison of the density of states near the band edge in 
the SCBA for the SRI and CI for u = 0.2eV and two val- 
ues of the impurity concentration. When ni = lO^^cm"^ 
we take V = 2eVnm^ and when rtj — lO^^cm^^ we take 
V' = l.GeVnm^. We see that the two scattering mecha- 
nisms with these parameters give almost identical results 
for the DOS. 



Appendix B: Optical conductivity 

In this appendix, we present the details of the calcula- 
tion of the optical conductivity in the presence of disor- 
der. The real part of conductivity at a frequency w can 
be calculated from the Kubo formula^ 



a{uj) = Re 



(Bl) 



where 11^^ (cj) is the Fourier transform of the real-time 
retarded current-current correlation function defined by 



1 1"°° 



Ja{t),MO) 



(B2) 



Here f2 is the area of the system and Jq is the electric 
current operator defined by 



Ja = (-e) 



(B3) 



where u^''^(k) = (kA', | Va |k, A), cj^^ and CkA are creation 
and annihilation operators for a state |A, k), respectively. 

A real-time retarded correlation function can be eas- 
ily calculated from a finite-temperature correlation func- 
tion through the analytic continuation in the complex 
w— space. The finite-temperature current-current corre- 
lation function is defined by 



I3hv 



A,A',k 

X 5kA(«l^n)3kA'(«'^n + '^'^n), (B5) 



where gkx{iujn) is the finite-temperature Green's func- 
tion defined by the Fourier transform of g]i\{T) — 

— (^TT-C]i\{T)cl.^{0)'^ . Note that the spectral representa- 
tion of g]i\{iujn) is given by 



5kA(iw„) 



duj' px{k,uj') 



27r iujr. 



■ UJ' 



(B6) 



where pA(k, w) is a spectral weight function with 
PA(k,a;) = — 2Im [5kA(w + i??)]. From Eqs. (|B5p and 
(lB6l). 



1, \ A / 



duj did'' 
V ^ I 2^"2^ 

k,A,A' 



.pA(k,c.')PA'(k,w") 



^ A'An \ AA'n \ /(^") ~ fi^') mv\ 



where f{uj) = [e^*^''"^^' + l] and fi is the chemical 
potential. Here the following frequency sums for fermions 
at finite temperatures were used: 



I3h 



1 ^ e*^"'' 

Ph ^ [iuJn — Uj'){iuin + iVn - tjj") 



(B8) 



Through the analytic continuation of w„ — > uj + irj, the 
conductivity expression can be derived from Eqs. (jBl[) 
and (|B7p . The energy bands have rotational symmetry 
and by taking this into account along with the spin and 
valley degeneracy, we finally get the conductivity expres- 
sion at zero-temperature as in Eq. ([9]). 



hllaisiil^n) 



fin 



dre^'^-^ {TrMr)MO)) , (B4) 
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